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ABSTRACT 

Context. AGB stars have long been held responsible for the important star-to-star variations in light elements observed in Galactic 
globular clusters. 

Aims. We analyse the main impacts of a first generation of rotating intermediate-mass stars on the chemical properties of second- 
generation globular cluster stars. The rotating models were computed without magnetic fields and without the effects of internal 
gravity waves. They account for the transports by meridional currents and turbulence. 

Methods. We computed the evolution of both standard and rotating stellar models with initial masses between 2.5 and 8 M Q within 
the metallicity range covered by Galactic globular clusters. 

Results. During central He-burning, rotational mixing transports fresh CO-rich material from the core towards the hydrogen-burning 
shell, leading to the production of primary I4 N. In stars more massive than M ^ 4M Q , the convective envelope reaches this reservoir 
during the second dredge-up episode, resulting in a large increase in the total C+N+O content at the stellar surface and in the stellar 
wind. The corresponding pollution depends on the initial metallicity. At low- and intermediate-metallicity (i.e., [Fe/H] lower than or 
equal to ~ -1.2), it is at odds with the constancy of C+N+O observed among globular cluster low-mass stars. 

Conclusions. With the given input physics, our models suggest that massive (i.e., > 4 M Q ) rotating AGB stars have not shaped the 
abundance patterns observed in low- and intermediate-metallicity globular clusters. Our non-rotating models, on the other hands, do 
not predict surface C+N+O enhancements, hence are in a better position as sources of the chemical anomalies in globular clusters 
showing the constancy of the C+N+O. However at the moment, there is no reason to think that intermediate mass stars were not 
rotating. On the contrary there is observational evidence that stars in clusters have higher rotational velocities than in the field. 

Key words, stars: AGB and post-AGB, rotation, abundances - globular clusters: general 
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' 1. Introduction Massive AGB stars that undergo efficient hot-bottom burn- 

J> , ing (HBB) during the thermal pulse phase (TP-AGB) have 

<3 ' While * e abundances of heay y elements (i-e-. Fe-group and a- been pr0 posed as the po ssible GC pol l uters in this so-called 

^ ■ elements) are fairly constant from star-to-star in any well-studied se lf-enrichment scenario (IVentura et al] l200lh . The AGB hy- 

H ; individual Galactic globular cluster (GCfl giant and turnoff stars pot hesis has been extensively discussed in the literature, 

, present some striking anomalies in their light elements content first on a qualitative basis, and more recently with the 

(Li, F, C, N, O, Na, Mg, and Al) that are not sh ared by their hel p of custom-made standard (i.e., non-rotating) stellar m od- 

field counterparts (for a review see |Gratton||2007|). Important is e ls (IVentura et alJllOOll |200H IVentura & p'Antonall2005 JEII3, 

that the sum [(C+N+0)/Fe] for individual GC stars otherwise 2008 albt iDenissenkov & Herwid |2003l iKarakas & Lattanziol 

presenting very different chemical features is constant to within |20Q3UHerwigl l2004allbt iFenner et alJl2004t iBekki et aO lofm 

experimental errors in all the clusters studied so far, with the These studies point out several difficulties in building the ob- 

possible exception of NGC 1851 (for references see §5). serV ed chemical patterns in theoretical TP-AGB models. The 

We now have compelling evidence that these peculiar chem- main problem stems from to the competition between the third 

ical patterns were present in the intracluster gas from which dredge-up (3DUP) that contaminates the AGB envelope with the 

second-generation (anomalous) low-mass stars formed, and that helium-burning ashes produced in the thermal pulse and HBB 

they resulted from the dilution of pristine material with the that modifies the envelope abundances via the CNO-cycle and 

hydrogen-burning products ejected by a first generation of the NeNa- and MgAl-chains. It is thus very difficult to obtain 

more massive and f aster evolving GC stars (for a review see simultaneous O depletion and Na enrichment in the TP-AGB 

iPrantzos et alj|200~7l) . envelope, while keeping the C+N+O sum constant as required 

by the observations (see Charbonnel 2007, for more details and 

Send offprint requests to: T. Decressin, email: decressin@astro.uni- references). 

bonn.de To date, only AGB models in the mass range 5-6.5 M G 

1 Except id Cen (Norr is & Da Cos ta 1995; Johnson et al. 20Q3). have managed to simultaneously achieve an encouraging agree- 
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ment with the observed O- depletion and the Na-enrichment 
dVentura&D'AntonalteOOSal). These models include the FST 
formulation for convection dCanuto et alJI 19961) . which strongly 
affects the O depletion on ce the stars enter the TP-AGB 
dVentura & D'Antonal l2005b). Compared to the classical MLT 
treatment, FST indeed leads to higher temperatures at the base 
of the convective envelope (resulting in more advanced nucle- 
osynthesis) and induces higher surface luminosities resulting in 
stronger mass loss and thus fewer thermal pulse s and 3DUP 
events. As shown bv lVentura & D'Antonal d2005bl) . some O de- 
pletion can also be obtained with the MLT prescription, but only 
when the free parameter a is arbitrarily and strongly modified 
with respect to the value calibrated on the Sun. These models 
are also able to produce a slight decrease in Mg accompanied 
with a large increase in Al. 

The Na-enrichment is more difficult to estimate and requires 
fine tuning of the NeNa-cycle reaction rates. More precisely, an 
increase in sodium is achieved in the most oxygen-poor ejecta 
of the 5-6.5 M models only when the maximum allowed values 
for the 22 Ne(p,y) rate are adopted. In s ummary, and as clearly 
stated by IVentura & D'Antona d2008al) . the AGB scenario is 
viable from the nucleosynthesis point of view, provided that 
only massive AGB stars of 5 to 6.5 M contribute to the GC 
self-enrichment, and under the physical assumptions described 
above. 

However the previously quoted models have only focused 
on physical uncertainties related to the TP-AGB phase, and 
their predictions have not been tested in different astrophysi- 
cal contexts. In particular, the impact of rotation in our mod- 
els on the nucleosynthesis predictions for AGB stars has never 
been investigated in the context of the GC self-enrichment sce- 
nario, although rotation is oft en invoked to understand a wide 
variety of observatio ns (e.g. iMaeder & Mevnetl l2000l l2006t 
IChiappini et alj 120081 see § 4). The present paper addresses 
this question for the first time, using up-to-date treatment for 
rotation-induced processes. As we shall see, the main signature 
of rotation on the chemical composition of the stellar envelope 
and winds of intermediate-mass stars already show up during 
the second dredge-up event and cannot be erased during the TP- 
AGB phase. 

2. Physical input of the stellar models 

Although only massive AGB stars in a very narrow mass range 
between 5 and 6.5 M Q are now suspected to play a r ole in 
the self-enrichment scenario dVentura & D , Antonal2008allb1) . we 
computed standard and rotating models of 2. 5, 3, 4, 5, 7, and 
8 M stars w ith the code STAREVOL (V2.75) dSiess et al.l2000t 
Siess 2006). We present results for several metallicities, namely 
Z = 4 x 10 \ 10~ 3 , 5 x 10~ 4 , 10~ 4 , and 10~ 5 (i.e., [Fe/H]^ -0.66, 
-1.26, -1.56, -2.26, and -3.26 respectivelyfl The com position is 
scaled solar according to the lGrevesse & Sauvalld 19981) mixture 
and enhancement in a-elements ([a/Fe] = +0.3 dex) is accounted 
for. All models were evolved up to the completion of the 2DUP . 

We used the OPAL opacity tables (Iglesias & Rogers! [T996I) 
above T > 8000 K that account for C and O enrichments, and 
the IPerguson et al. (2005]) data at lower temperatures. We fol- 



2 The models at [Fe /H]^ -0.66 and -3.26 are actually from 
Maeder & Meynej] d2001h and iMevnet & Maederl j2002h . respectively. 
To check the consistency of our predictions with those of the Geneva 
code, we computed a 7 M Q model at Z = 1CT 5 with STAREVOL and 
found excellent agreement in terms of both evolutionary and chemical 
characteristics (see Fig.O. 



lowed the evolution of 53 chemical spe cies from H to CI us- 
ing the NACRE nuc lear reaction rates dAngulo et al.lll999l) by 
default and those by | Caughlan & Fowled ([1988) otherwise (see 
ISiess & Arnouldl2008l) . The treatment of convection is based on 
the classical mixing length formalism with otmlt = 1.75, and no 
convective overshoot is i ncluded. The mass loss rates are com- 
puted with Reimers ( 1975) formula (with t;r = 0.5). 

For the treatment of rotation-induced processes we proceed 
as follows. Solid-body rotation is assumed on the ZAMS and 
a typical initial surface velocity of 300 km s _1 is chosen for 
all the models. (The impact of this choice is discussed in § 3.3 
where we present models computed with initial rotation veloc- 
ities ranging from 50 to 500 km s .) On the main sequence, 
the evolution of the internal angular momentum profile is ac- 
counted for with the comple te formalism developed by Zahn 
(119921) and IMaeder & Zahiil (1 19981) (see iPalacios et all I2003L 
2006; Decress in et al.ll2009l for a description of the implemen- 
tation in STAREVOL), which takes advection by meridional cir- 
culation and diffusion by shear turbulence into account. The ini- 
tial solid-body profile relaxes on the main sequence on a short 
timescale of a few Myr and generates differential rotation. For 
our most massive stellar models (M > 4 M©), the complete 
formalism for angular momentum transport is applied up to the 
2DUP. For the lower mass models, however, the complete treat- 
ment is only applied up to the end of the main sequence, while 
a more crude approach is used in the more advanced evolution 
phases where angular momentum evolves only through local 
conservation {i.e. only the structural changes modify the angular 
momentum). This simplification is motivated by the evolution- 
ary timescale (i.e., the Kelvin-Helmholtz timescale on the RGB) 
becoming shorter than the meridional circulation timescale. We 
test this assumption by running a 5 and 7 M Q with and with- 
out the full treatment of rotational mixing and the results show 
little difference, validating our approximation. In all cases the 
transport of chemical species resulting from meridional circu- 
lation and both vertica l and horizontal turbulen ce is computed 
as a diffusive process dChabover&Zahnl fl992) throughout the 
evolution. 

3. Signatures of rotation-induced mixing up to the 
early-AGB phase 

3.1. Models at [Fe/H]=-1.56 

Up to the beginning of the TP-AGB phase, the surface compo- 
sition of the standard models is only modified by the dredge-up 
event(s). Let us note that only the 2.5 and 3 M Q models undergo 
deep enough first dredge-up so as to modify their surface abun- 
dances; however, all stars experience to various extents the sec- 
ond dredge-up (hereafter 2DUP) after central He exhaustion. As 
described below, in the rotating models the effects of additional 
mixing become visible at the stellar surface during the quiescent 
central He-burning phase. 

Figure Q] shows the abundance profiles (in mass fraction) of 
'H, 4 He, 12 C, 14 N, and 16 0, as well as the total sum C+N+O, at 
the end of central He-burning and before the 2DUP, for standard 
(lower row) and rotating (upper row) 5 M model. 

In the standard case, 14 N steadily increases as one moves in- 
wards through the H-rich radiative layers down to the He buffer 
(located between M r ~ 0.7 and 1 .2 M Q ) as a result of CNO pro- 
cessing. Further inside, the stellar core has experienced complete 
He-burning and is essentially made of 16 and 12 C. In this 5 M Q 
model, during the 2DUP, the convective envelope penetrates into 
the He-rich buffer (down to M r ^ 0.85 M ) as indicated by 
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Fig. 1. Chemical profiles (in mass fraction) at the end of central He-burning in rotating (top) and standard (bottom) 5 M models at 
Z = 0.001, 0.0005, 0.0001, and 0.00001 (left to right). The elements shown are 4 He (dotted-dashed lines), 12 C (short-dashed lines), 
14 N (long-dashed lines), 16 (dotted lines), sum C+N+O (thick lines). The hatched area on top of each panel indicates the maximum 
extent of the convective envelope during the 2DUP. 



the hatched area. This produces an envelope enrichment in H- 
burning ashes: the surface abundances of 12 C and 16 decrease, 
while that of 14 N and 4 He increase. The sum C+N+O, however, 
remains constant, since only the H-burning products are dredged 
up. 

Rotation-induced mixing strongly modifies the internal 
chemical structure. As shown in the upper-row panels of Fig. [1] 
the abundance gradients are smoothed out in the radiative enve- 
lope (i.e., the region above the He-rich buffer up to the surface or 
the convective envelope) compared to the standard case: 14 N pro- 
duced in the internal H-burning layers diffuses outwards, while 
the 12 C and ls O present in the envelope are transported inwards. 
As a consequence, during the whole central He-burning phase, 
rotational mixing produces a continuous surface increase in 14 N 
concomitant to a decrease in 12 C and 16 0. At the same time, the 
products of central helium-burning, namely 12 C and 16 0, also 
diffuse outward in the He-rich buffer (i.e., the region where He 
is the dominant species, between 0.6 and 1.4 M in Fig.Q]). 

The transport of chemical species is mainly driven by shear 
turbulence (D S h ear ~ 10 8 cm s _1 ) in the radiative envelope and 
in the He-rich layer. At the interface of these two regions, i.e. at 
the base of the HBS, the large mean molecular weight gradient 
strongly reduces the efficiency of mixing (£> s hear ~ 10 3 cm s _1 ), 
thus preventing the transport of primary C, N, and O from the 



He-rich buffer to the HBS, hence thus to the surface during cen- 
tral He-burning. At the same time, hydrogen from the envelope is 
also transported inwards and rapidly captured by 12 C and 16 nu- 
clei through CNO burning at high temperature. This leads to the 
production of a peak of primary 14 N at the base of the hydrogen- 
burning shell (HBS) as seen in Fig. Q] The resulting chemical 
profiles at the end of central He-burning thus differ significantly 
from those obtained in the standard case where 14 N is only pro- 
duced in the HBS from the 12 C and ls O or iginally present in 
the st ar and is therefore of secondary origin ( Mevnet & Maeder 
2002). During the subsequent 2DUP, the convective envelope of 
the 5 M rotating model reaches the polluted He-rich buffer (see 
Fig. [TJ producing a large increase in 4 He. Simultaneously the 
primary CNO and thus the overall metallicity increases in the 
envelope. 

In the 2.5 and 3 M rotating models, the convective envelope 
does not reach the contaminated He-buffer during the 2DUP. As 
a consequence in these low-mass models, only the H-burning 
products are dredged up to the surface: 14 N strongly increases, 
while 12 C and 16 decrease, but the sum C+N+O, as well as 
the total metallicity, remain unchanged (see Table Q]). 4 He also 
diffuses outward into the radiative envelope leading to a surface 
4 He enhancement by about 0.03 (in mass fraction) compared to 
non-rotating models. 
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Table 1. Surface abundance variations after the completion of 
the 2DUP with respect to the initial composition (<5[X/Fe] = 
[X/Fe] 2DU p - [X/Fe] init ) for the models with initial value of 
[Fe/H]=-1.56 





2.5 M 


3M 


4M Q 


5M 


7 Mq 


8M Q 








Standard models 






He 


0.26 


0.26 


0.29 


0.32 


0.36 


0.36 


<5[C/Fe] 


-0.20 


-0.25 


-0.29 


-0.31 


-0.30 


-0.34 


5[N/Fe] 


0.44 


0.48 


0.60 


0.75 


0.79 


0.81 


c5[0/Fe] 


-0.01 


-0.01 


-0.03 


-0.06 


-0.10 


-0.11 


<S[CNO/Fe] 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 








Rotating 


models 






He 


0.29 


0.29 


0.32 


0.34 


0.35 


0.36 


<J[C/Fe] 


-1.12 


-0.86 


-0.20 


1.44 


2.14 


2.16 


<S[N/Fe] 


1.23 


1.12 


1.13 


1.24 


1.19 


1.09 


<5[0/Fe] 


-0.65 


-0.36 


-0.01 


0.18 


1.46 


1.55 


<S[CNO/Fe] 


0.00 


0.00 


0.15 


0.78 


1.62 


1.71 
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Fig. 2. Surface C+N+O increase index (<5 [C+N+0] 2 dup = 
[C+N+0/Fe] 2D up - [C+N+0/Fe] init ) at the end of the 2DUP 
for rotating stars with various initial masses and metallic - 
ities. Squares and triangles indicate models computed with 
STAREVOL and with the Geneva code, respectively (see text 
for details). 



Table[T]summarises the abundance variations after the 2DUP 
in all our standard and rotating models at [Fe/H] =-1.56. The 
main signature of rotational mixing at the surface of massive 
early- AGB stars (M > 4 M Q ) is a strong increase in He-burning 
products, i.e., primary CNO (see also Fig. |2). As all rotating 
models produce an enrichment of CNO elements inside the He- 
buffer, the total increase of C+N+O at the surface of rotating 
models mainly depends on the depth reached by the convective 
envelope during the 2DUR We thus obtain a stronger variation 
with increasing stellar mass. 




200 400 

< V surf> ( km S_1 ) 



600 



Fig. 3. Top: enrichment in C+N+O after the 2DUP completion 
as a function of the mean main sequence velocity for a 5 M Q , Z = 
0.0005 for initial surface velocities of 0, 50, 150, 300, 500 km s 1 
as indicated. Bottom: dis tribution of surface velocities observed 
bv lMartavan et all (120061) in the young LMC cluster NGC 2004 
for B-type stars with mass from 2 to 10 M Q (v sin i measurements 
are multiplied by a factor 4 /if). 



3.2. Influence of metallicity 

The impact of rotation on stellar properties and stellar 
yields is known to depe nd strongly on metallicity (see, e.g., 
iMevnet & Maeden |2002|) . The metallicity dependence of our 
predictions is depicted in Fig. [TJ showing the internal profiles of 
12 C, 14 N, and le O at the end of central He-burning in standard and 
rotating 5 M models at four different metallicities ([Fe/H] 
-1.26, -1.56, -2.26, and -3.26). As previously explained, the 
C+N+O profile outside the CO core is constant in the standard 
models, while it strongly increases in the He-rich layers below 
the HBS in the rotating models. This C+N+O step is higher in 
lower metallicity stars, which results in a stronger CNO surface 
enrichment after 2DUP in the most metal-poor stars as shown 
in Fig. [2] The maximal depth reached by the convective enve- 
lope for a given stellar mass hardly depend on metallicity; e.g., it 
reaches 1 .04 and 1 .07 M for the 7 M at Z = 10 3 and Z = 10 -4 . 

At [Fe/H] =-2.26, the envelope (and thus the wind) of all the 
stars more massive than ~ 5M Q undergo a C+N+O increase by 
1 to 2 orders of magnitude, while an increase by a factor of 5 is 
obtained at [Fe/H] = -1.26. When metallicity becomes higher 
than [Fe/H] -1, rotation-induced mixing increases the total 
C+N+O by less than a factor 2-3, which is undetectable with 
current observations (see § 5). The effect is null at [Fe/H] = 
-0.66. 



3.3. Influence of initial rotation velocity 

We present several models for the 5 M Q star at Z = 0.0005 
([Fe/H] = -1.56) with initial rotation velocities between and 
500 km s _1 as part of investigating the dependence of our results 
on this parameter. As shown in Fig. [3] the C+N+O enhance- 
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ment increases with increasing initial rotation. A higher initial 
velocity leads to a faster spinning core and to a larger differential 
rotation at the core edge during central He-burning. Since mix- 
ing stems mainly fro m shear turbulence and scales as (dCl/dr) 2 
dTalon&Zahdl 19971) . more CNO elements are stored in the He- 
rich buffer and then revealed at the surface after 2DUP. In the 
model with an initial velocity of 150 km s -1 , the C+N+O rises 
by a factor 3, while it increases by more than a factor 14 for an 
initial velocity of 500 km s _1 . 



4. Discussion on the model uncertainties 

4.1. Mean rotation velocity 

As shown above, the total C+N+O enhancement depends on the 
assumed initial rotational velocity. We note that, for an initial 
(i.e., ZAMS) velocity of 300 km s , the time-averaged velocity 
of our models on the main sequence ranges between 220 and 
256 km s depending on the stellar mass and metallicity. What 
are the observational constraints on this ingredient? 

iMartavan etall d2006l l2007al) find mean vsinz' of 161 ± 
20 km s" 1 and 155 + 20 km s' 1 for SMC B-type stars of 2-5 M 
(111 stars) and 5-10 M (81 stars), respectively. Accounting in 
a statistical way for the sin / effect (i.e. multiplying the aver- 
aged vsinz by 4/n supposing a random distribution of the ro- 
tational axis), we obtain averaged values for v between 197 and 
205 km s for B-type stars in the SMC. At first glance, the main- 
sequence time-averaged velocity of our models is slightly higher 
than the observed values for SMC stars. However, the following 
points have to be considered: 

1 . The mean observed values quoted above do not account for 
Be-type stars. Let us recall that Be stars are fast rotators that 
present emission lines originating in an outward equatorial 
expan ding disk probably for med due to strong stellar rota- 
tion (IMartavan et al.ll2007bl) . Therefore Be stars belong to 
the upper part of the velocity distribution, and it is legitimate 
to incorporate them in the estimate of the averaged veloci- 
ties of B-type stars. Taking them into account does enhance 
the observed average velocities. As an illustrative example, 
the mean v sin i for SMC Be stars in the mass ranges 2-5, 5- 
10 M are 277 ± 34 km s" 1 (14 stars) and 297 ± 25 km s" 1 (8 1 
stars), respectively. As can be seen in Fig. [3] our assumptions 
on the stellar ro t ation v elocities ar e thus totally realistic ; 

2. IMartavan etall (l2007al) (see also iHunter et al.ll2008bl) find 
that, for both B and Be stars, the lower the metallicity, the 
hi gher the rotational ve loci ties. This agrees with the finding 
by lMaeder etall d 19991) and lWisniewski & Bjorkmanl d2006h 
that the fraction of Be stars with respect to the total num- 
ber of B+Be stars increases when the metallicity decreases. 
Since the metallicities considered in this work are lower than 
that of the SMC, we may expect that the averaged velocities 
of the stars would be somewhat higher than the one quoted 
above for the SMC; 

3. Last but not least, there is evidence that the rotation rates 
of stars are higher in clusters than in the field (lKellerll2004t 
Huang & Giesl 120051; IStrom et al J 120051: iDufton et alj 120061 
Wolff et al.ll2007l) . 

Thus, in view of these remarks, it does appear that our adopted 
rotation velocities are probably very close to the averaged veloc- 
ities of stars in clusters for the range of masses and metallicities 
considered. If the previous extrapolations are correct, the values 
given in Fig. [2] should be considered as lower limits. 



4.2. Treatment of rotation 

Of course, rotating models are not free from uncertainties 
and their predictions should be carefully compared with well- 
observed features. The physics included in the present models 
is the same as the one that provides a good fit to the following 
observed features: 

- the s urface enrichments in nitrogen in main-sequence B-type 
stars (Maeder et al. 2009), even if invoking binarity or mag- 
netic fi elds is required to e xplain the whole observational 
pattern (IHunter et al.ll2008al) : 

- the ob served number ratio of b lue to red supergiants in the 
SMC dMaeder & Mevnetll200l : 

- the variation with t he metallicity of the number fr action of 
WR to O-type stars dMevnet & Maederj|2003l , l2005l) : 

- the variation with metallicity of the nu mber ratio of type 
Ibc to type II core collapse supernovae dMevnet & Maederl 
120051) . 

- the lithium abundance patterns in A-type and early F- 
type dwarf stars, as well as in their subgiant descendants 
dCharbonnel & Talonlll999HPalacios et alj|2003l) . 

Moreover, they provide a reasonable explanation for the ori- 
gin of the high N/O observed at t he surface of metal-poor hal o 
stars, as well as for the C/O upturn (Chiappini et alj2006l[2008b . 
Therefore, while the present models are not free of uncertainties, 
they have the nice property of accounting for the above observed 
features. 

4.3. Rotational rate of remnants 

The present models have some difficulty, however, in account- 
ing for the obser ved rotat i on rates of young pulsar s and white 
dwarf s (see e.g. iKawalerl 119881; iHeger et all 120051; ISuiis et al.1 
2008). More precisely, they predict too fast rotation of the stellar 
cores in the advanced phases. This may be stem from different 
causes that could lead to additional angular momentum loss from 
the central regions at different evolutionary phases: 

- already during the nuclear lifetime; 

- at the time of the supernova explosion in the case of neutron 
stars or during the TP- AGB phase at the time of the super- 
wind episode in the case of white dwarfs; 

- by the neutron stars or the white dwarfs themselves shortly 
after their formation. 

In the second and third cases mentioned above, the physics and 
the predictions of the present models would not need to be re- 
vised since the loss of angular momentum would occur after 
the evolutionary stages covered by the present computations. 
However, if the loss of angular momentum occurs before sec- 
ond dredge-up, another mechanism should be included at earlier 
phases in rotating models. 

The study of the s-process nucleosynthesis in low-mass TP- 
AGB stars could provide some hints to the rotational evolution of 
the stellar core. At the moment the properties and the behaviour 
of rotating TP-AGB stars are poorly known. Calculations us- 
ing adiffusiv^tteatment for the transpo r t of angular momentum 
bv lLanger et al.1 dl999l) . iHerwig et al] d2003h . and ISiess et all 
(2004) have shown that rotationally induced instabilities provide 
enough mixing to trigger the 3DUPs. However, a shear layer at 
the base of the convective envelope leads to an efficient pollu- 
tion of the 13 C pocke{3 (the neutron source) by 14 N with the re- 

3 In this framework, the neutrons needed for the s-process are re- 
leased in the l3 C pocket by C(or,n) reaction. 
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suit of strongly inhibiting the s-process nucleosynthesis. This is 
at odds with the observations of s-stars and tends to indicate that 
the modeling of rotation used in these models needs to be im- 
proved and that angular momentum must have been removed by 
the time the thermal pulses start. Therefore, if the extraction of 
angular momentum occurs before core helium burning, no CNO 
enrichment will be expected because it relies on shear due to a 
fast spinning core. On the other hand, if the extraction occurs 
during the early-AGB phase, it could allow both the CNO en- 
richment of the He-buffer and the s-process during the TP- AGB 
phase. 

The mechanism frequently pointed to for removing angu- 
lar mo mentum from the core is the magnetic field. For in- 
stance, iHeger et ail d2005l) have shown that magnetic coupling 
between the core and en velope can account for the rotation rates 
of young neutron stars. ISuiis et aTl ((2008) also find that mag- 
netic torques may be required to understand th e slow spin rate 
of white dwarfs. However the dynamo model (jSpruit 2003) on 
which the current rotating models with magnetic fields are based 
ha s recently b e en te sted through hydrodynamical computations 
by IZahn et all (120071) . who do not find the amplification of the 
magnetic field as expected from the theory, casting some doubt 
on its validity. In view of the remaining theoretical uncertainties 
associated with the treatment of magnetic fields, it seems reason- 
able to stick to the present models whose predictions account for 
a broad variety of observations as described in § 4.2. 

Finally let us note that iTalon & Charbonnel d2008l) predict 
that angular momentum transport by internal gravity waves 
should be efficient in intermediate-mass stars during the early- 
AGB phase. If these waves were the culprit, then the predictions 
of the present models would be valid as far as the CNO enrich- 
ment is concerned since it builds up earlier in the life of the star. 



5. Consequences for the self-enrichment scenario 

5. 1 . Summary of the theoretical predictions 

During central He-burning, rotational mixing efficiently trans- 
ports primary 12 C and 16 outside the convective core in the 
H-burning region where these elements are processed by the 
CNO cycle, resulting in an important production of primary 
14 N. In the massive (M <: 4M ) rotating models, after central 
He-exhaustion, the convective envelope penetrates into the lay- 
ers affected by rotation-induced mixing. In contrast to standard 
models, the 2DUP produces a large surface enrichment in total 
C+N+O that cannot be erased by hot bottom (hydrogen) burn- 
ing during the subsequent TP-AGB evolution. Third dredge-up 
episodes further increase the total C+N+O mass fraction as they 
bring the products of He-burning to the surface. 

As a consequence, if rotating massive AGB stars were re- 
sponsible for the abundance patterns observed in GCs, one 
would expect large C+N+O differences between (O-rich and Na- 
poor) first-generation star s and (O-poor and Na-rich) s econd- 
generation stars (see e.g. iPrantzos & Charbonnelll2006l) . Such 
differences would easily amount to 1 .6 dex in GCs with [Fe/H] — 
1.5, and to 2.2 dex in the most metal-poor clusters with [Fe/H] — 
2.2 for the cases where the ejecta of the first-generation would 
not have been diluted with pristine material. Diluting these ejecta 
with 10 times more pristine material would still give differences 
of 0.6 to 1.2 dex, well above the observed dispersion (see § 5.2). 
Actually, to maintain C+N+O constant would require such a 
high dilution of the ejecta with the intracluster matter that the 
O-Na and Mg-Al anticorrelations would be erased and, as far 



as abundances are concerned, first and second generation stars 
would be indistinguishable. 



5.2. Comparison with the observations 

C, N, and O abundances have been determined simultaneously 
in stars of several GCs. Up to now, no significant star-to-star 
variation of the total C+N+O has been detected, except in NGC 
1851. 

In the case of t he most metal-ric h GCs such as 47 Tuc 
(ICarretta et al.1 12001 and NGC 6712 (lYong et al.l l2008h . this 
C+N+O constancy is actually compatible with pollution by both 
standard and rotating AGB stars since the 2DUP CNO enrich- 
ment is found to be negligible in the most metal-rich rotating 
models. This is of course without considering possible C+N+O 
increase due to the 3DUP of He-processed material during the 
TP-AGB phase, which is beyond the scope of the present study. 

In the case of intermediate-metal licity ([Fe/H1 ~ - 1-2) 
GCs such as NGC 28 8 and NGC 362 Jbickens et alj [l99ll) and 
M4 divans et all 1 19991: ISmith et alj |2005h. the total C+N+O is 
found constant from star to star, within observational uncertain- 
ties that amount to +0.25-0.3 dex. At this metallicity, the most 
massive intermediate-mass rotating models (M > 5) predict a 
C+N+O increase between 0.5 and 0.7 dex at odds with the ob- 
servations. The only exception is NGC 1851 (Yong et al. 2009, 

[Fe/H] 1.2), where variations in the total CNO by 0.57 ± 

0. 15 dex were found, which is compatible with the CNO enrich- 
ment produced by our rotating models as shown in Fig. [2] We 
note, however, that the stars observed by Yong and collaborators 
in NGC 1851 are very bright objects. In view of its position in 
the colour-magnitude diagram, the only star with C+N+O ex- 
ceeding the typical error bar could actually be an AGB star. In 
that case, the high value of the total CNO may not have been 
inherited at birth by the star, but may instead be due to nuclear 
processes within the star itself. Being an AGB would also ex- 
plain why this star is the only one among the sample to exhibit 
some enhancement in s-process elements. We note that in this 
cluster a double sub giant branch has b een detected by photomet- 
ric measurements (Milone_etal][2008) that can be interpreted as 
caused by an age difference am ong cluster stars of about 1 GyrQ. 
However. ICassisi et al.l (|2008) propose to fit the double sequence 
assuming two coeval cluster stellar populations with a C+N+O 
difference of a factor of 2. CNO measurements in NGC 1851 
unevolved stars (as done in most of the other clusters studied so 
far) are thus crucial for settling the problem. 

Finally, the CNO predictions for the rotating most mas- 
sive intermediate mass models lead to increasing factors 
up to 10 (5 M ) and 100 (7 M Q ) at low-metallicity 
([Fe/H] = -2.2), in clear contradiction with the observations 
in the metal-poor clusters studied so far, i .e., M3 and M13 
dSmith et al .1119961: ICohen & MelendezlfeoOSl) . NGC 6752, and 
NGC 6397 dCarrettaetalJl2005lP t CNO measurements in other 
very metal-poor globular clusters like M15 or M92 would be ex- 
tremely valuable in this context. 



4 The hypothesis of two populations with similar age and [Fe/H] dif- 
fering beyond the err orbars have been r uled out by the spectroscopic 
study of RGB st ars bv|Yong et al. (2009) and by the photometric study 
of RR Lyrae bv lWalkej dl998l) . 

5 NGC 6397 ([Fe/H] 1.95) does not show a very extended O- 

Na anticorrelation, a more modest C+N+O increase is expected in that 
case. 
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6. Conclusions 

In this paper we have investigated the effects of rotation in low- 
metallicity intermediate-mass stars during their evolution up to 
the completion of the 2DUP and arrival on the AGB. Our rotat- 
ing st e llar m odels include the complete fo rmalism developed by 
Zahn (119921) and iMaeder & Zahnl d 19981) and accounts for the 
transport of chemicals and angular momentum by meridional 
circulation and shear turbulence. With respect to standard mod- 
els, the most important change concerning the nucleosynthesis is 
the large 14 N production resulting from the diffusion of protons 
below the HBS in a region enriched with primary C and O during 
central He-burning. During the subsequent 2DUP, the convective 
envelope of massive AGB stars deepens in this region, producing 
a large surface increase in the total C+N+O, which irreversibly 
imprints the yields. 

This behaviour is in sharp contrast to what is observed in 
low- and intermediate-metallicity globular clusters where the 
sum C+N+O is constant with in the expe rimental errors. Our ro- 
tating models based on the IZahnl (I 19921) formalism, which ne- 
glect the effects of magnetic fields and internal gravity waves, 
thus suggesting that massive rotating AGB stars can be dis- 
carded as potential polluters in the self-enrichment scenario in 
globular clusters, unless the crowded environment prevented 
intermediate-mass stars from rotating. This latest hypothesis is 
highly improbable in view of the observations finding higher 
stellar rotation velocities in clusters than in the field. 
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